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Abstract 



It has recently been observed that MgB2 is a superconductor with a high 
transition temperature. Here we propose a model of anisotropic s-wave su- 
perconductivity which consistently describes the observed properties of this 
compound, including the thermodynamic and optical response in sintered 
MgB2 wires. We also determine the shape of the quasiparticle density of 
states and the anisotropy of the upper critical field and the superfluid density 
which should be detectable once single-crystal samples become available. 
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1. Introduction 

The recently discovered superconductor MgB2 has received great attention, mainly be- 
cause of its high transition temperature of approximately 39K. This relatively high Tc 
value may not be too surprising since MgB2 involves some of the lightest elements in the 
periodic table, thus implying a large Debye phonon frequency for this compound. Indeed, 
data from thermodynamic measurements suggest a Debye frequency ©d ~ 750 K as 
well cLS cL SI gnificant isotope effect with a = 0.26. 

More surprisingly, a rather small zero-temperature energy gap Aq = 1.20kBTc was de- 
duced from measurements of the specific heat and from the optical conductivity [^. 
Presently, we do not know of any other examples for such a small ratio Aq/ ksT^ in conven- 
tional s-wave superconductors. This observation naturally suggests an anisotropic s-wave 
order parameter for this material, where the energy gap detected in the thermodynamic 
measurements is in fact the minimum of the superconducting gap function. An alternative 
explanation has been suggested in terms of multi-band models. In many respects the 
proposed phenomenological model with an anisotropic energy gap is similar to a multi-band 
model. In fact, the latter may be thought of as a discretized version of the former. More 
recently, the upper critical critical field for a dense wire of MgB2 was inferred from the field 
dependence of the critical current. P| From these measurements it was concluded that MgB2 
is an extreme type II superconductor with a Ginzburg-Landau parameter k ~ 23. 

Drafting an effective model with an anisotropic s-wave order parameter for MgB2, the 
first important question to address is the direction and the magnitude of the anisotropy 
over the approximately ellipsoidal Fermi surface. In the distantly related tetragonal systems 
YNi2B2C and LuNi2B2C, such an anisotropy has been shown to appear in the a-b planes. 



T0|] On the other hand, MgB2 has a hexagonal crystal structure, and thus the anisotropy 
is most likely to occur along the c-direction similar to other hexagonal crystals, such as UPts 
Tl| and UPd2Al3 [|1^], and tetragonal crystals, such as Sr2Ru04 |T3|JI^ . 



We therefore propose a BCS model for MgB2 with a superconducting order parameter 
given by 



Afk) 
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1 + az^ 
1 + a 



(1) 



where the parameter a determines the anisotropy, z = cos 6, and 6 is the polar angle. For 
the calculations outlined in this paper, we impose the experimentally observed gap ratio, 

Amin{T = 0) Ao 



1.20, 



(2) 



fcijTc (1 + a)kBTc 

which in turn yields a ~ 1. In the following, we will therefore explore the consequences of 
anisotropic s-wave superconductivity in MgB2 within the framework of weak-coupling BCS 
theory, setting a = 1. 
2. Density of States 

In Fig. 1(a) the anisotropic s-wave order parameter A(k) = A(l + z^)/2, is plotted in 
momentum space. This function is an ellipsoid with a minor axis of length A/2 within the 
a-b plane, and a major axis of length A along the c-direction. The smaller magnitude of the 
gap function within the a-b plane is consistent with the stronger in-plane Coulomb repulsion 



that has been suggested in Ref. 115 1. 
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FIG. 1. (a) Anisotropic s-wave order parameter, (b) Quasiparticle density of states. 



The corresponding density of states can be calculated from A(k) within weak-coupling 
theory. It is given by 



N{E) 



dz^e 



\E\ 



E2- A2(l + z2)2/4^ 
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if < E < A/2, 
if A/2 < ^ < A, 

if E > A, 



(3) 



where and A;) are the complete and incomplete elliptic integrals of the first kind. 

This quasiparticle density of states is shown in Fig. 1(b). It is fully gapped with an onset 
of spectral weight at the minimum gap value E = A/2 and a cusp aX E = A. 

3. Thermodynamics 
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FIG. 2. (a) Energy gap in an anisotropic s-wave superconductor, (b) Entropy (solid line) and 
specific heat (dashed line). Above Tc these quantities are equal (dot-dashed line), (c) Thermody- 
namic critical field for anisotropic (solid line) and isotropic (dashed line) s-wave superconductors. 
Inset: deviation of these critical fields from a parabolic temperature dependence. 



Following the formalism by Bardeen, Cooper, and Schrieffer [0, we calculate the tem- 
perature dependence of the energy gap A(T), the entropy 5'(T), the specific heat C(T), and 
the thermodynamic critical field Hc{T) for the anisotropic s-wave superconductor. Here, the 
entropy is given by 

S{T) = -A dEiV(E)(/ln/ + (l-/)ln(l-/)), (4) 

•J 
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(5) 



where / = (1 + expf3E) ^. Furthermore, the specific heat and the upper critical field are 
obtained from 



C(T) = T 



dS{T) 
dT 



and 



in 



dTS{T). 



(6) 



These results are displayed in Fig. 2. As expected, the characteristic jump of the specific 
heat at Tc, AC/Cn = 1.18, is small compared with the value 1.43 for an isotropic s- 
wave superconductor. This is consistent with the experimental picture for MgB2. 
Furthermore, it is seen in Fig. 2(b) that the thermodynamic critical field is also reduced 
in the anisotropic case. The deviation of the critical field from a parabolic temperature 
dependence, defined by 



D{T/T,) = 



Hc{T) 



(7) 



is shown in the inset of Fig. 2(c). The magnitude of this deviation is substantially larger 
than for the isotropic case, again consistent with the experiments. H] 




FIG. 3. Superfluid densities of anisotropic and of isotropic s-wave superconductors. 

Finally, the superfiuid density Ps(T) shown in Fig. 3 is found to be rather anisotropic. 
Similar to the energy gap and the thermodynamic critical field, upon increasing the tem- 



perature its departure from its zero-temperature value is exponentially small in the low- 
temperature regime. On the other hand, in the vicinity of the superfluid density for the 
isotropic case is linear in temperature, whereas for the anisotropic case Ps{T) is observed to 
be slightly concave along the c-direction and slightly convex in the a-b plane. Once single 
crystals of MgB2 become available, this anisotropy in the superfluid density should serve as 
a signature for the order parameter symmetry. 
4- Upper Critical Field 

Recent measurements of the upper critical field in sintered samples of MgB2 have reported 
a somewhat unusual temperature dependence. This is likely related to the apparent 
anisotropy of the superconducting order parameter, and should thus become even more 
evident in measurements of the upper critical field in single crystals. Therefore we consider 
here the two cases of an external magnetic field parallel and perpendicular to the crystal 
c-axis of MgB2. 
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FIG. 4. Upper critical fields of anisotropic and isotropic s-wave superconductors. 



(a) H II c 



The upper critical field along the crystal c-direction is determined from 0,0] 



(8) 
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where p = {y%eHc2{T)) j {2{2iiTY) and z = cos 9. In Fig. 4, the solution of this integral 
equation normalized by its derivative at hc2iT) = Hc2(T) / {—dHc2(T) / dT)]^^, is plotted 
along with the /ic2(7')-curves for the a-b-direction and for the isotropic case. In the limit 
T ^ we obtain /ic2(0) — 0.824 which is somewhat larger than the corresponding value for 
the isotropic s-wave superconductor, /ic2(0) — 0.728. In the opposite limit, T h'^2iT) 
exhibits a rather linear temperature dependence within the weak-coupling BCS model. This 
behavior is also observed for the isotropic s-wave superconductor. 
(b) H II a 

For the upper critical field within the a-b plane a mixing occurs of the zeroth and the 



second Landau levels, leading to [p!0|,p!7 



A(k,r)~A(l + C(at)2)|0), 



(9) 



where |0) is the Abrikosov state, is the raising operator, and C is the mixing coefficient 
between the Landau levels which has to be determined self-consistently along with the critical 
field. This leads to a set of coupled integral equations, 
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exp 



{-pu\i-z^))y (11) 



where again p = {v'^peHc2{T)) / {2{2-nTY), and sit? = 1 — z^. The numerical solution of 
this set of equations is shown in Fig. 4. The mixing coefficient (not shown) is found to 
decrease monotonously from 0.069 down to 0.062 as the temperature is increased. This 
implies a relevant admixture of the second Landau level. In the zero-temperature limit, 
we obtain h'^2i^) — 0.664. By comparing /ic2(0) with /ic2(0), we can thus conclude from the 
model calculation that the critical field observed in sintered MgB2 represents H^2{T) because 
H^2iT) > H^2iT) for all temperatures. We also note that H^2iT), H^2iT) > Hc{T) since for 
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sintered MgB2 wires the Ginzburg-Landau parameter was observed to be very large k ^ 23. 
i 

5. Conclusions 

Based on presently available experimental data on MgB2 we have constructed a model 
of anisotropic s-wave superconductivity for this compound with Amm/^max = 1/2- This 
simple theory appears to account rather well for the thermodynamic properties |^,|], the 
optical measurements, and for the upper critical field data on sintered samples of 

MgB2. Experiments on single crystal samples will clearly be important to further investigate 
the applicability of this model. 
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